Pseudo-spin-dependent scattering in carbon nanotubes 
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The breaking of symmetry is the ground on which many physical phenomena are explained. This 
is important in particular for bipartite lattice structure as graphene and carbon nanotubes, where 
particle-hole and pseudo-spin are relevant symmetries. Here we investigate the role played by the 
defect-induced breaking of these symmetries in the electronic scattering properties of armchair single- 
walled carbon nanotubes. From Fourier transform of the local density of states we show that the 
active electron scattering channels depend on the conservation of the pseudo-spin. Further, we show 
that the lack of particle-hole symmetry is responsible for the pseudo-spin selection rules observed 
in several experiments. This symmetry breaking arises from the lattice reconstruction appearing at 
defect sites. Our analysis gives an intuitive way to understand the scattering properties of carbon 
nanotubes, and can be employed for newly interpret several experiments on this subject. Further, 
it can be used to design devices such as pseudo-spin filter by opportune defect engineering. 



I. INTRODUCTION 

The unique electronic properties of single-walled car- 
bon nanotubes (SWNTs) — related to their unusual 
band structures 1,2 — have attracted great attention in 
fundamental and applied research because of the possi- 
bility of exploring phenomena unique to one-dimensional 
systems. 3,4 Defectless SWNTs are well characterized 
both theoretically and experimentally — however, the 
presence of defects or tube endings is in fact crucial 
for the observation of quantum mechanical coherence 
phenomena. 5 7 In this respect, scanning tunneling mi- 
croscopy /spectroscopy (STM/STS) represents a powerful 
instrument of investigations. This is an unparalleled tool 
able to measure the local electronic properties of SWNTs 
and correlate these with the atomic structure of the tube. 
Thanks to this technique, it has been possible to visual- 
ize Friedel oscillations 8 10 and quasi-bound-states 11,12 in 
metallic SWNTs. These articles, reconstructing the low- 
energy scattering spectrum of SWNTs, have routinely 
shown the presence of an asymmetry in the scattering 
properties of left- and right-moving electrons. However, 
only in the work by Ouyang et al 9 there is an attempt 
to explain this asymmetry by relating it to the nature of 
defects. Their general conclusion is that this asymme- 
try roots in the different symmetry properties of the 7r 
and 7r* bands in the armchair SWNT compared to the 
symmetry of defects. 

Here, we present an in depth study of the scattering 
properties of electrons in armchair SWNTs in presence of 
structural defects. We show how the scattering proper- 
ties are altered when fundamental symmetries of SWNTs, 
as the pseudo-spin symmetry and particle-hole symme- 
try (PHS), are broken by the defects. In particular we 
focus on di-vacancy (DV) and Stone- Wales (SW) defects 
(c.f. Fig. 1). These are prevalently observed in graphene 
and SWNTs, 2,3,13 16 however our approach is quite gen- 
eral and can be extended to other types of defects. We 
show that pseudo-spin symmetry is determining the ac- 
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FIG. 1: The DV defect in Panel (a) (Type I) is symmetric 
with respect to an exchange of the two carbon atoms A and 
B in each unit cell of the honeycomb lattice, whereas for the 
DV defect in Panel (b) (Type II) this symmetry is broken. 
In Panel (a), the atoms 3 and 4 are missing in the atomic 
structure while the bonds 1-2 and 5-6 represent the lattice 
reconstruction. Panel (c) is a SW defect. (See main text for 
further details.) 

tive scattering channels. Further, contrary to the hy- 
pothesis by Ouyang et al 9 , we substantiate that the ex- 
perimentally observed asymmetry of electron scattering 
arises from the breaking of PHS caused by the recon- 
struction of the chemical bond between carbon atoms at 
the defect site. PHS breaking is also at the origin of other 
experimental observations 10 12 unexplained so far. This 
newly revealed feature of defects in bipartite lattices can 
also be exploited to design a defect-based device able to 
filter the pseudo-spin species. 



II. MODEL AND FORMALISMS 

We consider an armchair SWNT in presence of DV 
and SW defects. Simulations 14 16 and experiments 13 
indicate that DVs are likely the prevailing defect-type 
in SWNTs and graphene. Additionally, DV defects in 
SWNTs energetically favor the so-called 585 reconstruc- 



2 




-x r x 



k 

FIG. 2: Complete energy spectrum of an armchair SWNT 
[Panel (a)] and its linearization around the charge neutrality 
point [Panel (b)]. In Panel (b), the full arrows represent In- 
traVB scattering, dashed arrows are InterVF scattering and 
dotted arrows are Inter VB scattering. The dashed and solid 
lines represent electronic states with different pseudo-spin a. 

t ion 1 5,1 6 where two new chemical bonds forms, therefore 
creating two pentagons and an octagon [c.f. Figs. 1(a) 
and 1(b) — named here Type I and Type II, respectively]. 
SW defects consist of a local tt/2 rotation of a C-C bond 
[c.f. Fig. 1(c)], creating two pentagons and heptagons. 
We describe the SWNT as a finite stripe of graphene 
with periodic boundary conditions along the boundaries 
delimited by the chiral vector of the tube. We treat the 
7r-orbital electrons within a tight-binding approximation. 
The defectless Hamiltonian with spin-independent hop- 
ping term t reads %q = — 1 5^<i j> c J c j? wnere the sum 
runs over all pairs of nearest-neighbor carbon atoms and 
c\ [ci] is the creation [annihilation] operator for a p z elec- 
tron at atom i. Two different methods are used in order 
to study armchair SWNTs: an approximate low-energy 
one where only the linear dispersing modes are considered 
— named Method A (c. f. App. A), and a exact numerical 
method within tight-binding model — named Method B. 



Here, the SWNT is placed seamlessly between two semi- 
infinite leads at the same chemical potential consisting as 
well of SWNTs. Within Method A, we evaluate modifica- 
tions to the scattering properties induced by defects via 
Fermi's Golden Rule 18 (c. f. App. A). With Method B, 
the local density of states (LDOS) V(x, E) is evaluated 
as the imaginary part of the retarded Green's functions 20 
(c. f. App. B). Then for each energy value E a Fourier 
transform (FT) is performed on the LDOS for the coordi- 
nate along the tube axis in order to obtain the FT-LDOS 
D(2fc, E). Here, the factor 2 originates from the fact that 
|^(x)| 2 is probed. 

Within the tight-binding approximation the Hamilto- 
nian Hdv describing a 585 DV is obtained by removing 
the hopping terms to the missing carbon atoms and by 
introducing hopping elements to for the newly formed 
carbon bonds. It is important to note that these new 
bonds link carbon atoms of the same sub-lattice. For the 
case of SW defects, two distinct reconstruction hopping 
elements tj^ 1 and t^ 2 (tj} 1 < t < £d 2 ) are introduced into 
the Hamiltonian Hsw- 

In general a defect can give rise to three different types 
of scattering processes [c. f. Fig. 2(b)]: inter- valley for- 
ward/backward (Inter VF/B) and intra- valley backward 
(IntraVB). However, only Inter VB processes are pseudo- 
spin conserving. 

In the following we consider the two possible orienta- 
tions of a 585 DV defect on an armchair SWNT shown 
in Fig. l(a)-(b), Type I and II, respectively, and the SW 
defect in Fig. 1(c). 



A. Type I di- vacancy defects 

The Hamiltonian associated with this defect com- 
mutes with the pseudo-spin symmetry operator S 
(c. f. App. CI), therefore the pseudo-spin a is a good 
quantum number. Within Method A, the scattering ma- 
trix element of the perturbation Hbv between different 
SWNT eigenstates {\ctFk)} reads 



(aFK,\HBv\<r'F'K') = -^-U B cos[(A + + 5+)R x ](l + aa') + 2t cos 
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Here A± = F ± F f , 5± = k ± and R x = x/3a C c/2. 
The conservation of the pseudo-spin a is also immediately 
evident from Eq. (1) since only pseudo-spin conserving 
scattering processes (a = cr') have a non vanishing am- 
plitude. Therefore, only Inter VB processes are induces 
by Type I DV defects. Equation (1) yields two different 
energy dependent scattering probabilities V{E) for states 



with pseudo-spins a = ± =t / 4- 



a^a(E) - 



, 2tt ER x 

t d — 2crt cos h cr 

3 7 



(2) 



Here E = ±7^ with 7 = 3acct/2 is the energy disper- 
sion around the charge neutrality point (CNP). Without 
lattice reconstruction (£d = 0) the two scattering proba- 
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FIG. 3: Unnormalized scattering probabilities as a function 
of the energy for different scattering events. Panel (a) refers 
to Type I DV defect and Panels (b) and (c) to Type II DV 
defect. In all Panels red and blue lines refers to pseudo-spin 
conserving processes of I and t carriers, respectively. In Panel 
(a) we show only InterVB processes, where the dashed and 
the solid lines refer to £d = and £d = 0.5 1, respectively. In 
Panels (b) and (c) the solid lines refer to InterVB scattering, 
the dashed-black lines to InterVF scattering and the dotted- 
black lines to IntraVB scattering. The lattice reconstruction 
parameter is £d = in Panel (b) and £d = 0.5 1 in Panel (c). 
The value of tr> is obtained from DFT calculations 19 . In Panel 
(d) we display the case of the SW defect with tr> 1 = 1.2t and 
t Dl = o.St. 



bilities are particle-hole symmetric 

However, the lattice reconstruction breaks PHS 
(c. f. App. C2). For £ D ^ and energy values around 
the CNP, we observe that the InterVB scattering process 
1<=>1 is strongly suppressed compared to the pro- 
cess. This is illustrated in Fig. 3(a) where the full and 
the dashed lines refer to t& 7^ and £d = 0, respectively. 
Hence, we can clearly associate the different scattering 
probabilities to the breaking of PHS associated with the 
lattice reconstruction. 

Using Method B, we consider a (5,5) armchair SWNT 
with a single Type I DV defect. In Fig. 4 we display the 
FT-LDOS V(2k, E) for the case without [panel (a)] and 
with lattice reconstruction t^> = 0.5t [panel (b)]. For en- 
ergies within the range ±1.5 eV we find a direct image of 
the linear SWNT spectrum around E = 0. This signa- 
ture corresponds to InterVB scattering as explained by 
the extended k • p model 11 ' 12 . In the absence of lattice 
reconstruction the FT-LDOS is particle-hole symmetric 
[Fig. 4(a)], 

V(2k,E) = V(2k,-E) 

whereas the lattice reconstruction introduces a strongly 
reduced scattering probability for the a =1 scattering 
channel [c.f. Fig. 4(b)] compared to the a =t scattering 
channel around the CNP. 



B. Type II di-vancancy defects 

The Hamiltonian associated with the tilted Type II 
defect does not commute with the pseudo-spin symme- 
try operator S in contrast to the Type I defect. In this 
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FIG. 4: Density plot of the Fourier transformed local den- 
sity of states as a function of the exchanged momentum and 
energy for the case of a single Type I, Type II DV and a SW 
defect. In Panel (a) we have a Type I DV defect (£d = 0) 
whereas in Panel (b) £d = 0.5t. In Panel (c) we have a Type 
II DV defect with t D = 0.5 t. In Panel (d) we have a SW 
defect with £d x = 1.2£ and tr> 1 = O.St. Note that for energies 
higher than ±|1.5| eV the nonlinear SWNT-subbands appear. 



case we obtain more intricate results for the scattering 
processes. Method A shows that the pseudo-spin non- 
conserving InterVF and IntraVB processes have a finite 
probability [c.f. Figs. 3(b) and 3(c)]. Again, without lat- 
tice reconstruction the scattering is particle-hole sym- 
metric, i.e., 

V^(E)=V^(-E) 
V a ^- a (E)=V a ^- (J (-E) 

hold [c.f. Fig. 3(b)]. However for t B ^ PHS is bro- 
ken again — we observe, contrary to the previous case, 
a strong enhancement of the scattering probability for 
the ^=>^ InterVB processes whereas the InterVB 
scattering is almost suppressed [c.f. Fig 3(c)]. This is 
confirmed also by Method B. In Fig. 4(c) we show the 
FT-LDOS for the case of a Type II DV defect with a 
finite value for the lattice reconstruction (tr> = O.bt). 
Compared to the case of Type I DV defect are now clear 
the IntraVB scattering processes at k > nm" 1 . 



C. St one- Wales defects 

The SW defects intrinsically include a lattice recon- 
struction — therefore PHS is broken by its very nature. 
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For an infinite SWNT there is a degeneracy for the SW 
defect related to clockwise and ant i- clockwise rotations 
of the C-C bond. In order to correctly describe these de- 
fects, we introduce a symmetrized Hamiltonian account- 
ing for both configurations. This Hamiltonian commutes 
with the pseudo-spin symmetry operator <S, therefore also 
in this case the pseudo-spin <j is a good quantum num- 
ber. Indeed, we find by using Method A that also for the 
SW Hamiltonian only scattering processes conserving the 
pseudo-spin are allowed [Fig. 3(d)]. 

In Fig. 4(d) we present the case of a single SW defect 
in a (5,5) SWNT within Method B. Here again there are 
no signs of processes where pseudo-spin is not conserved. 
However, opposite to the DV cases, the V(2k, E) shows 
a more complex dependance on the energy and pseudo- 
spin. 



III. CASE OF TWO DI- VACANCIES 

We investigate — only within Method B — a (10,10) 
armchair SWNT with two reconstructed Type I and II 
DV defects (£d = 0.5t) placed at a distance of L = 
14.5 nm along the SWNT axis. Since t D ^ 0, PHS is al- 
ready broken. Four unequal combinations can be formed 
from the two Type I and II defects — however we con- 
sider here only two relevant cases as depicted in the insets 
of Fig. 5. The scattering behavior of each of the single 
defects determines the pattern of V(2k, E) between the 
defects. The overall features of the FT-LDOS for the 
combination of two Type I DV defects [Fig. 5(a)] resem- 
bles those of a single Type I DV defect shown in Fig. 4(b). 
We observe, however, a discretization in energy indicat- 
ing the formation of quasi-bound-states between the de- 
fects. The pseudo-spin-selection rule for Type I DV de- 
fects accounts for the fact that those bound states are 
essentially only formed by a =t states as the asymme- 
try in the amplitude between the different linear slopes 
clearly shows. Figure 5(b) contains one Type II DV 
defect that breaks the pseudo-spin symmetry — now In- 
terVF and IntraVB scattering occur. The fingerprint of 
the Inter VF processes are the vertical lines at k ~ 9 nm -1 
and ~ 17 nm -1 . The IntraVB scattering is revealed by 
the dispersive features at k > nm -1 . 



IV. THEORY VS. EXPERIMENTS 

We now apply these findings to interpret the experi- 
mental results of Ouyang et a/. 8,9 and Lee et al. 10 . In 
the first work the scattering at a single isolated defect 
in an armchair SWNT is analyzed. The reconstructed 
energy spectrum shows a selection rule of Inter VB scat- 
tering processes with respect to the pseudo-spin. In view 
of our results, we can attribute the asymmetry to the 
lattice reconstruction. In the second case the standing 
waves of electrons scattered at one end of a nanotube 
peapod are examined. Also in this case a clear asymme- 




FIG. 5: Density plot of the Fourier transformed local den- 
sity of states as a function of the exchanged momentum and 
energy for the case of two 585 DV defects with tD = 0.5t. In 
Panel (a) we consider two Type I DV defects , in Panel (b) two 
Type I and Type II DV defects. The signal at k ~ 25.5 nm" 1 
is associated with R x , the smallest length scale of the arm- 
chair SWNT along the x-axis. 



try in the scattering processes is found. We can relate 
this asymmetry to the peapod at the tube end. There 
various hexagons of the lattice structure are replaced by 
pentagons therefore mixing the two sub-lattices of the 
SWNT as in the case of the DV and SW defects. Thus 
the breaking of PHS explains the asymmetry observed 
experimentally. In the experiments by Buchs et al. 11 the 
quasi-bound-states in defected metallic SWNT were in- 
vestigated. There the defects were constituted mainly 
by DVs created by means of ions irradiation 12 . In order 
to make a comparison we have analyzed the case of two 
DV defects in two possible spatial configurations. Also 
in this case we have observed an absence of processes 
that are not conserving the pseudo-spin in the case of 
two defects Type I. For the second configurations where 
processes not conserving pseudo-spin are allowed, the nu- 
merical outcome accounts quite well for the experimental 
results. 

Summarizing these finding we recognized two main fea- 
tures affecting the electron scattering at defects in the 
lattice structure of an armchair SWNT: the conservation 
of the pseudo-spin degree-of-freedom, and the breaking 
of PHS due to the lattice reconstruction. The former 
determines the active scattering channels, the latter in- 
troduces an asymmetry among the different scattering 
processes. 



V. PERSPECTIVE 

The methods we have introduced here can also be em- 
ployed to study defects in graphene and chiral SWNTs. 
In the latter case, the pseudo-spin is in general not con- 
served. However, the breaking of PHS due to the lat- 
tice reconstruction is a general feature of the underlying 
SWNT honeycomb lattice and does not depend on the 
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tube chirality. 

Choosing the symmetry of the defect accordingly, it 
is possible to modify the ratio of the transmission prob- 
abilities of the two pseudo-spin species. Therefore, the 
prospect of engineering definite defects in SWNTs, would 
permit to filter electrons with a specific degree-of-freedom 
- the pseudo-spin — in analogy to a similar effect on 
valley-spin induced by line defects in graphene 21 . 

As we have shown with our results, the pseudo-spin de- 
pendent scattering response changes considerably with 
the electron energy. Moreover, different defects can be 
combined in a line, so as we have demonstrated, in order 
to obtain an enhanced effect of filtering. Therefore, by 
an opportune defect sequencing and by fixing the appro- 
priate energy it should be possible to select the desired 
pseudo-spin in a controlled way. The problem of engi- 
neering well defined defects, e.g. via ion irradiation, can 
be overcome by employing adatoms instead of DV and 
SW defects. These can be placed on the SWNTs with 
an atomic precision using, e.g. the tip of a STM. In fact, 
adatoms break PHS in the same manner as the DV and 
SW defects by allowing hopping between tt electrons of 
carbon atoms of the same sub-lattice 22 . 



Appendix A: Method A 

In proximity of the charge neutrality point the energy 
spectrum of armchair SWNTs can be approximated by 
two linearly dispersing branches centered around the two 
valleys K$ = ±47r/3-\/3acCo where acc ~ 0.142 nm is 
the C-C bond length (c.f. Fig. 2). In this approximation 
the electronic wave function can be expressed as 17 



(r\aFK) 



9=A,B 



fp*<P P F K (r) . 



(Al) 



Here F = zbifo is the valley index, n is the momentum 
relative to F, and a is the pseudo-spin. This is defined by 
a = sgn(6F), where b = ± = Right/Left is the electron 
motion direction. The coefficients f pa account for the 
two inequivalent carbon atoms (p = A,B) in the lattice 
structure. These are defined as 



J_fl P = A 
y/2 \cr p = B ' 



(A2) 



The function (/} p ]? K (r) describes the armchair SWNT 
wave function component on sub-lattice p. The functions 
^PpFk{t) in Eq. (Al) is defined as 



1 



(F+K)R a 



X(r-r p ), 



where the sum runs over the SWNT lattice sites, 2Nl 
is the total number of carbon atoms, and x( r ~ r p) are 
the localized p z orbitals at the positions r p = Ft + t p on 
sub-lattice p, and t p the displacement vector within the 



honeycomb lattice unit cell 17 . Finally, R x is the projec- 
tion of the lattice vector onto the tube axis (c.f. Fig. la). 

Within Fermi's Golden Rule approximation, the tran- 
sition probability between different states {\<jFk)} is 
given by 

Vf^i (x \ f (aFK\H dsEe ct\aFK) i \ 2 S(E f -E i ) 

where E^/f are the energies of the initial and final states, 
respectively. 



Appendix B: Method B 

The local density of states (LDOS) has been evaluated 
numerically via the retarded Green's functions 



0(E) 



1 



E - H + S R + S L 



(Bl) 



Here T~L is the Hamiltonian of the defected SWNT, and 
E R / L are the self-energies of the right or left lead, re- 
spectively. They have been determined by a recursive 
Green's function method. 20 From G(E) we have evalu- 
ated the LDOS T>(x, E) at carbon atom site x via 



V(x,E) 



1 



lm[(x\g(E)\x}} 



(B2) 



where Im[. . .] is the imaginary part. In order to get a 
smooth LDOS on the surface of the SWNT, we have per- 
formed a convolution of V(x, E) with a function of the 
form exp(— Ar) where r is the distance of a fictitious tip of 
a STM from the carbon atom at x, and A is an opportune 
constant. 



Appendix C: The symmetries 

1. The pseudo-spin symmetry 

In our case the pseudo-spin operator S is given by the 
Pauli matrix a x . In general, electron scattering is pseudo- 
spin conserving if the defect Hamiltonian commutes with 
S. Since the effect of S is only a simultaneous exchange 
of the two carbon atoms of the unit cell, see Fig. 1(a), 
pseudo-spin is conserved if and only if the defect geom- 
etry is symmetric with respect to an exchange of the A 
and B atoms within the single unit cells. Hence it is 
clear that the Hamiltonian associated with the Type I 
DV defect commutes with <S, whereas this is not true for 
the case of the Type II DV defect. In the case of the 
symmetrized version of the SW defect in Fig. 1(c), its 
Hamiltonian commutes with the pseudo-spin operator. 

2. The particle-hole symmetry 

The particle-hole symmetry operator for the SWNT 
is defined as IA = — o z . A Hamiltonian T~L is particle- 
hole symmetric if the following properties hold: \p[h]) = 
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U\h[p)) where is a particle [hole] state, and 

{T-LM} — 0- When the lattice reconstruction is bond- 
ing together carbon atoms of the same sub-lattices, the 
Hamiltonian associated with the reconstructed defect is 
not anti-commuting with U. 
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